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Abstract 

' We carefully analyse the present experimental evidence for a com- 

plex CKM matrix, even allowing for New Physics contributions to 
^ , ex, aj/q,x s , AMs d , AMb 3 , and the AI = 1/2 piece of B — ► pp and 

| B —> pir. We emphasize the crucial role played by the angle 7 in 

both providing irrefutable evidence for a complex CKM matrix and 
placing constraints on the size of NP contributions. It is shown that 
CN \ even if one allows for New Physics a real CKM matrix is excluded at a 

99.92% C.L., and the probability for the phase 7 to be in the interval 
J§ ! [-170°; -10°] U [10°; 170°] is 99.7%. 

1 Introduction 

At present, the Cabibbo-Kobayashi-Maskawa (CKM) P mechanism for 
>• ■ flavour mixing and CP violation is in agreement with all available experimen- 

^ . tal data. This is a remarkable success, since it is achieved with a relatively 

small number of parameters. Once the experimental values of |V^ S | , \ V c b\ and 
\V u b\ are used to fix the angles Qyi-.Q^z and 6*13 of the standard parametriza- 
tion, one has to fit, with a single parameter 5 13 , the experimental values of 
a large number of quantities, including e K , sin (2/9) , AM Bd , as well as the 
bound on AMb s . This impressive result is nicely represented in the usual 
unitarity triangle fits |2j. In view of the remarkable success of the Standard 
Model, it is plausible that the CKM mechanism gives the dominant contri- 
bution to mixing and CP violation at low energies, although there is still 
significant room for New Physics (NP). 

1 On leave of absence from Departamento de Fisica and Centro de Fisica Teorica de 
Particulas (CFTP), Instituto Superior Tecnico, Av. Rovisco Pais, 1049-001 Lisboa, 
Portugal. 
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In this paper, we address the question of how the data available at present 
already provide an irrefutable proof that the CKM matrix is "non-trivially 
complex", thus implying that the charged weak-current interactions violate 
CP. Obviously, in the framework of the Standard Model (SM), the CKM 
matrix has to be complex, in order to account for the observed CP violation, 
both in the kaon and in the B sector. The above question only becomes non- 
trivial if one allows for the presence of NP j^J. We shall carefully analyse 
the present experimental indications in favour of a complex CKM matrix, 
correcting at the same time various misleading and in some cases erroneous 
statements one finds in the literature. 

To have an irrefutable proof that the CKM matrix is complex is of the 
utmost importance in order to investigate the origin of CP violation, as well 
as to analyse what classes of theories of CP violation are viable in view of 
the present experimental data provided by B-factories. For example, one of 
the crucial questions concerning the origin of CP violation is whether CP 
is violated explicitly in the Lagrangian through the introduction of complex 
Yukawa couplings or, on the contrary, CP is a good symmetry of the La- 
grangian, only spontaneously broken by the vacuum. There are two classes 
of theories with spontaneous CP violation: 

i) Those where the phases arising from the vacuum lead to a complex 
CKM matrix, in spite of having real Yukawa couplings. In this class 
of theories, there are in general more than one source of CP violation, 
namely the usual KM mechanism, together with some new sources of 
CP violation. 

ii) Models where all phases can be removed from the CKM matrix and 
thus CP violation arises exclusively from New Physics. 

Examples of class i) are, for instance, two Higgs doublet extensions of 
the SM with spontaneous CP violation, without natural flavour conservation 
(NFC) in the Higgs sector |H E] , as well as models where CP is broken at 
a high energy scale jH] through a phase in the vacuum expectation value 
(vev) of a complex singlet, which in turn generates a non-trivially complex 
CKM matrix through the mixing of isosinglet quarks with standard quarks. 
Examples of class ii) are the simplest supersymmetric extensions of the SM 
with spontaneous CP violation [Zj, as well as three Higgs doublet models 
with spontaneous CP violation [8J and NFC in the Higgs sector. 

It is clear that if it can be proved that experimental data constrain the 
CKM matrix to be complex, even allowing for the presence of New Physics, 
the only models of spontaneous CP violation that are viable belong to class 
i). 
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The paper is organized as follows. In the next section we specify what 
our assumptions about NP are and we address the question of how to obtain 
evidence for a complex CKM matrix which would not be affected by the pres- 
ence of this type of NP. In section 3 we discuss the possibility of obtaining an 
experimental proof that the CKM matrix is complex, from the measurement 
of four indepedent moduli of this matrix. We argue that the cleanest proof 
would be obtained if one would use only moduli from the first two rows of 
the CKM matrix. However, we point out that contrary to some statements 
in the literature, this proof cannot be obtained in practice, since it would 
require totally unrealistic precision in the measurement of \V US \, \Vd\- In 
section 4 we analyse the impact of the measurement of AM Bd and clj/^k s 
and introduce a convenient parametrization of NP. In section 5 we emphasize 
the importance of a measurement of 7 either direct or through the combined 
measurement of (3, a. Finally, we present our summary and conclusions in 
section 6. 



2 Evidence for a complex CKM matrix unaf- 
fected by the presence of New Physics 

Our assumptions about NP will be the following. We will assume that in 
weak processes where the SM contributes at tree level, NP contributions are 
negligible, while NP may be relevant to weak processes where the SM only 
contributes at the loop level. This is a rather mild assumption on NP. In 
addition, for the moment we will also assume that the 3x3 CKM matrix 
is unitary, but we will come back to this point later. An example of this 
framework is the minimal supersymmetric standard model (MSSM). 

This scenario implies that NP could be competing with any SM loop and 
more likely in K°-K°, B G d -B% B° s -B° s and D°-D° mixing. Therefore, the 
usual interpretation of ex, clj/^>k S i AMs d , AMb 3 in terms of 613 gets invali- 
dated. The same could also be the case for the direct CP- violating parameters 
e'/e, A^(B° K + tt~) 9\ and various semileptonic rare decays. To clarify 
which are the fundamental observables that will not be polluted by NP con- 
tributions, let us write the CKM matrix in a particular phase convention ^0] 
in terms of all the rephasing invariant quantities fTJ E] that can a priori be 
measured in the flavour sector: 

/ \V ud \ \V us \e^' \V ub \e-^\ 
Vckm=[ -\V cd \ \V CS \ \V cb \ (1) 

V \Vtd\e~^ -\V u \jx \V tb \ ) 
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where the CP-violating phases introduced in Eq. are defined by: 



p = a r g (-v cd v; b v; d v tb ) 
x = Mg(-v ts v t tv*v; b ) 



i = ^g(-v ud v: b v; d v cb ) , 
x' = ^ g (-v cd v c * s v: d v us ) . 



(2) 



Note that a = arg (—V td V t * b V* d V ub ) obeys the relation a = ir — (3 — 7, by 
definition. Without imposing the constraints of unitarity, the four rephasing 
invariant phases, together with the nine moduli are all the independent phys- 
ical quantities contained in Vckm- In the SM, where unitarity holds, these 
quantities are related by a series of exact relations which provide a stringent 
test of the SM [HI [Hi- 
It is clear that the moduli of the first two rows are extracted from weak 
processes at tree level in the SM. From top decays we have some direct infor- 
mation on \Vtb\ - useless in practice compared to its unitarity determination 
from \Vub\ and |V^,| -, and essentially no direct information on \ V ts \ and \ V t d\ 
from tree level decays. Therefore, the extraction of the moduli of the third 
row is made with loop processes and consequently can be affected by the 
presence of NP. From Eq. it is evident that the phases (3 and x wm 
only enter in loop amplitudes because they appear in transitions involving 
the top quark, and therefore the usual extraction of these phases could be 
also contaminated by NP. x' is too small to be considered (see references 
llOp. Therefore, 7 is the only phase that can be measured without NP 
contamination, since it enters in b — > u transitions not necessarily involving 
loop mediated processes. 

In our framework, the most straightforward irrefutable proof of a complex 
CKM could come from the knowledge of the moduli of the first two rows of the 
CKM matrix and/or from a determination of 7 in a weak tree level process. 
To use the information on /3, x and the moduli of the third row will require a 
much more involved analysis including the presence of NP. But by the same 
token, this analysis will be extremely interesting in the not less important 
goal of discovering the presence of NP. In our search for an irrefutable proof 
of complex CKM, we will also analyse in parallel its consequences for the 
detection of the presence of NP ^3] , confirming the recent results by Ligeti 



3 Obtaining Im (V^V^V^V^) = J from four 
independent moduli 



Within this class of models, in order to investigate whether the present exper- 
imental data already implies that CKM is complex, one has to check whether 
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any of the unitarity triangles is constrained by data to be non-"flat", i.e. to 
have a non-vanishing area. If any one of the triangles does not collapse to 
a line, no other triangle will collapse, due to the remarkable property that 
all the unitarity triangles have the same area. This property simply follows 
from unitarity of the 3x3 CKM matrix. The universal area of the unitar- 
ity triangles gives a measurement of the strength of CP violation mediated 
by a ^-interaction and can be obtained from four independent moduli of 
Vckm- The fact that one can infer about CP violation from the knowledge 
of CP-conserving quantities should not come as a surprise [TK] . It just re- 
flects the fact that the strength of CP violation is given by the imaginary 
part of a rephasing invariant quartet [T7], J = ±Im {V ia VjpV*pV? a ) , with 
(i 7^ j, a 7^ /3), which in turn can be expressed in terms of moduli, thanks to 
3x3 unitarity. Restricting ourselves to the first two rows of Vckm, to avoid 
any contamination from NP, a possible choice of independent moduli would 
be \V US \, \V cb \, \V u b\ and \V c d\. One can then use unitarity of the first two 
rows to evaluate J, which is given, in terms of the input moduli, by 

4J 2 = A(l-\V ub \ 2 -\V us \ 2 )\V ub \ 2 \V cd \ 2 \V cb \ 2 - 

- (IK/ - \V cd \ 2 + \V cd \ 2 \V ub \ 2 - \V cb \ 2 \V ub \ 2 - \V cb \ 2 \V US \ 2 ) 2 . (3) 

Note that Eq. (JHJ) is exact, but the actual extraction of J from the chosen 
input moduli, although possible in principle, it is not feasible in "practice" . 

To illustrate this point, let us consider the present experimental values 
of \V US \ , \ V cb \ , \ V ub \ and \V cd \, assuming Gaussian probability density distri- 
butions around the central values. We plot in Fig. ^the probability density 
distribution of J 2 , generated using a toy Monte Carlo calculation [THj.Onlv 
31.1% of the generated points satisfy the trivial normalization constraints 
and, among those, only 7.9% satisfy the condition that the unitarity trian- 
gles close (J 2 > 0). 

In order to extract information on J 2 from the input moduli, it may be 
tempting to plot only the points that satisfy all unitarity constraints, i.e. 
normalization of columns and rows of Vckm ; together with the constraint of 
having J 2 > 0. We have plotted these points in Fig. |21 which naively lead 
to J = (2.6 ±0.8) x 1(T 5 . 

We wish to emphasize that, contrary to what has been stated in the 
literature ^H], one cannot interpret the above result as providing evidence 
for a complex CKM matrix, using as input data only |V^ S | , \ V cb \ , \ V ub \ and 
\V cd \. This can be trivially seen from the standard CKM parametrization 
[201121], by noting that from the input moduli \V US \ , \ V cb \ and \V ub \, Si 2 ,s 2 3 
and S13 can be obtained. In order to conclude that \V cd \ has necessarily a 
non- vanishing contribution proportional to cos 8\z (the leading one sensitive 
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(a) The complete distribution. 



(b) The region J 2 ~ 0. 



Figure 1: 
0.15)l(r 2 ; 



J 2 distribution from \V US \ = 0.2200 ± 0.0026, \V cb \ = (4.13 ± 
\V ub \ = (3.67±0.47)10- 3 and \V cd \ = 0.224 ± 0.012. 
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Figure 2: J 2 and J distributions. 



to S\s), one would need to know \V c d\ with a relative error of < 10~ 4 ! It 
is interesting to point out that the leading contribution to \V c d\ comes from 
\V US \, therefore also \V US \ would have to be known with the same level of 
precision, but not \V U b\ and \V c b\. For completeness we plot in Fig. |3(a)| the 
value of J 2 extracted from input moduli that are assumed to have the required 
unrealistic precision. It is clear that now, essentially all points have J 2 > 0, 
and therefore the corresponding J = (3.2 ± 0.5) x 10~ 5 has full meaning. 
Figure 3(b)| shows the J 2 distribution obtained with the same inputs as in 



Fig. |3(a) except for the central value of \V c d\; in this case only 0.16% of the 
points have J 2 > 0. Notice that both values, \V cd \ = 0.21985 ± 0.00002 and 
\V c d\ = 0.21955 ± 0.00002, are fully compatible with present measurements 
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(\V c d\ = 0.224 ±0.012). These two figures clearly show the required precision 
in \V US \ and \V cd \ to check 3x3 unitarity and to have an irrefutable proof of 
a complex CKM. 



j 1 



5. • 10 10 1. ■ 10 ' 1.5 • 10 ' 2. • 10 ' 



(a) \Vcd\ = 0.21985 ±0.00002 



j 1 
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(b) \V cd \ = 0.21955 ±0.00002 
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Figure 3: J 2 distributions assuming the unrealistic precision values \V US \ = 
0.22000 ± 0.00002 and \V cd \. 

Finally, it should be also emphasized that although the above analysis 
was made using the almost collapsed unitarity triangle corresponding to the 
first two rows of Vckm, the same conclusions can be obtained if we consider 
the "standard triangle", corresponding to orthogonality of the first and third 
columns. One can evaluate any of the angles of the unitarity triangle in terms 
of the input moduli, obtaining for example: 

2 7 \ ' * 11 \ ' I ' rl '\ I ' rlt \ I ' l ' ; / i\ 

sin — = ■ — ■ n n i , 4 

2 41^11^11^11^1 1 ^ 

Since we are using as input moduli | V^, s | , |Vrf>|, \V ub \ and \V c d\, it has to be 
understood that in Eq. (JH) one has: 

\v tb \ 2 = 1 - \v cb \ 2 - \v ub \ 2 , 

\V td \ 2 = \V us \ 2 -\V cd \ 2 + \V ub \ 2 , (5) 
\V u d\ = 1 — \V US \ — \V ub \ . 

The resulting expression can be easily written in terms of the input moduli 

as: 

\V cd \ 2 (1 - \V ub \ 2 ) - |Ks| 2 (1 - \V cb \ 2 ) + \V ub \ 2 \V cb \ 2 
cos 7 = -. . (6) 

2 \V ub \ \Vcdl \V cb \ \/l — \V US \ — \V ub \ 
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It is obvious that a "legitimate" extraction of 7 from the input moduli would 
require an unrealistic precision in the determination of the chosen moduli. 

I 2 1 2 

In particular, the difference \V c d\ — \V US \ would have to be known at the 

1 1 2 1 1 2 a 

level of \V US \ \V c b\ , thus requiring \ V c d\ and \V US \ to be measured with a 10 
relative error. This confirms our previous argument, based on the standard 
parametrization. Of course, Eq. (0J) can be used to determine sin 2 7/2 if 
one uses as input moduli \V US \, \V u t,\ ,\V c b\ and \V t d\- However, \V t d\ has the 
disadvantage that its extraction from experimental data is affected by the 
possible presence of NP. We will turn to this question in the next section. 

Above, we have argued that, although possible in theory, it is not viable 
in practice to prove that Vckm is non-trivially complex, from the knowledge 
of four independent moduli belonging to the first two rows of Vckm- We 
emphasize that this would be the ideal proof, since the extraction from ex- 
periment of the moduli of the first two rows is essentially immune to the 
possible presence of NP. 

4 The impact of AAfg d an d dj/^K s an d the 
NP parametrization 

The experimental measurement of the B^-B® oscillation frequency has reached 
an average accuracy of almost 1% . The theoretical prediction of AMg d 
within the SM is given by: 

G 2 M 2 

AM Bd = 2 |M 12 | = -JL^ VBd m Bd B B J 2 Bd S (x t ) \V tb V td \ 2 , (7) 

where we have followed the standard notation. Equation (J7J) is an approxi- 
mation valid up to order S (x c ,x t )/S (x t ) ~ 1CT 3 . From this expression we 
can extract |V^ 6 V^|. Furthermore, unitarity and the experimental values of 
\V u b\ and \Vcb\ determine \V t b\, thus allowing the extraction of |K<f|- In Fig. |U 
we plot the new J 2 distribution obtained with \ V t d\, thus determined together 
with \V US \, \Vcb\ and \V u b\. Table El presents all the values used for the input 
parameters. 

There is a striking difference between Fig. 0] and Fig. ^ As we have 
emphasized in the previous section, the majority of points obtained for J 2 , 
when one uses as input \V US \, \V c b\, \V u b\, \V c d\, do not satisfy the unitarity 
condition that the uc triangle closes (i.e. J 2 > 0). On the contrary, when we 
use \V US \, \V c b\, \V u b\, \V t d\ to extract J 2 , the majority of the points satisfies 
the unitarity condition that J 2 > and in fact the value of \J\ agrees with its 
experimental value \J\ ~ 1CT 5 (of course, using exclusively moduli of Vckm 
one cannot determine the sign of J). This means that, within the SM, one can 
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Figure 4: J 2 distribution obtained with \V US \, \V c b\, \V u b\ and \ V tc i\ (extracted 
from AM Bd ). 

use information on \ V u b \ / \ V c b\ and AM^ to conclude that the triangle cannot 
collapse to a line, thus implying a complex CKM matrix. Unfortunately, this 
is not an irrefutable proof that CKM is complex, since AM^ can receive 
significant contributions from NP, as we have emphasized in section El The 
same is true for the CP asymmetry aj/^ Ks in the B%, B® — > J/^K S decays. 
Let us introduce the following parametrization of the NP contribution to 
B^-Bj mixing: 

M 12 = r 2 d e- l2 ^[M 12 } SM (8) 

where [M 12 ]sm is the SM box diagram contribution. The NP contribution to 
M 12 implies that clj/^k s no longer measures 2/5, but one has instead: 

oj/*x s = sin 2(f3 - (j) d ) = sin 2(3 . (9) 

Although the expression for [M 12 ]sm is the one given by the SM, its actual 
numerical value may differ from the SM prediction, since models beyond the 
SM allow in general for a different range of the CKM matrix elements. 

It is clear that allowing for the presence of NP in Mi 2 , parametrized as 
in Eq. (JHJ), one can fit the data on aj/^^ s , AMb 4 , even with a real CKM 
matrix, by putting: 

(3 = 0; sin(20 d ) = -aj^ Ks (10) 

and then adjusting to fit AM^. Obviously, two solutions are obtained 
for corresponding to the two ways the unitarity triangle can collapse (i.e. 
7 = or 7 = 7r). It must be stressed that to write Eq. one has to assume 
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that the potential NP contamination to the penguin diagram contributing to 
the decay amplitude is negligible. This process is dominated by a tree and a 
penguin with the same phase; therefore, in this particular case, it is natural 
to assume no NP in the decay amplitude. For very small deviations from 
these assumptions see 

Independently of the question of whether the CKM matrix is complex or 
not, it is interesting to investigate what the allowed values for r d , cf) d are, 
taking into account the present data, summarized in Table El We calculate, 
using again the Monte Carlo method, probability density distributions for 
CKM parameters, r\ and 2<p d using as constraints Gaussian distributions for 
moduli of the first two rows of the CKM matrix, AM^ d and aj/q,K s - In 
Fig. |5(a)| we plot 68% (black), 90% (dark grey) and 95% (grey) probability 
regions of the probability density function (PDF) of the apex — V ud V* b /V cd V* b 
of the db unitarity triangle. In Fig. |5(b)| we represent joint PDF regions in 
the plane (r^, 2<p d ). Because 7 gives the apex of the triangle, it is clear 
from Fig. |5(a)| that there is essentially no restriction on 7. On the contrary, 
because the moduli of the first two rows put an upper bound on |/3| and 
upper and lower bounds on R t = \V td V t * b \ / \ V cd V* b \, we can see in Fig. |5(b)| 
significant constraints on 2<p d and r\. Although the experimental value of 
the semileptonic asymmetry Asl (the asymmetry in the number of equal 
sign lepton pairs arising from the semileptonic decay of B^-B® pairs) has 
some impact on these figures, we have not included it in order to show in a 
clear way the impact of the actual measurements. We will come back to this 
point later on. 

At this stage, the following comment is in order. We are assuming a class 
of theories beyond the SM that give NP contributions to AM^, but keep 
Vckm unitary. The most important example of this class of theories are the 
supersymmetric extensions of the SM. Since in this framework Vckm is still 
unitary, once \V US \, \V c b\, \V u b\ / \ V c b\ are fixed by data, there is only (3 as a 
free independent CKM parameter. It is then clear that with only two exper- 
imental inputs, namely AM# d , aj/q,K s , one cannot fix the three parameters 
r di 4>d, ft- One may be tempted to include information on AM^ s . However, 
in our scenario, this would require the introduction of new parameters (r s , 
<j) s ) giving NP contributing to AMg s . 

In the next section, we will emphasize the importance of 7 in obtaining 
an irrefutable proof that Vckm is complex, even allowing for the presence of 
NP. 
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Figure 5: 68%(black), 90%(dark grey) and 95%(grey) probability regions. 



5 The importance of a direct measurement of 
the phase 7 

In our framework, we can distinguish two conceptually different ways of mea- 
suring 7 in B® non-leptonic decays: those that come from the interference 
of two tree-level decays such as b — > sue with b — > scu, or b — > <iwc with 

6 — > dcit, and those other processes where the presence of both tree and pen- 
guin diagrams is allowed, such as b — > duu interfering with b — > d. The first 
group includes the usual methods that are called, in the literature, meth- 
ods to measure 7 or 2/3 + 7 |2lH 12*^] . In general they involve the transitions 
B — > Di^ or B ^ £)7r as the flavour-representative examples of the con- 
sidered decays. The second group includes the usual methods to measure a 
|2"H| EE] , or more generally a = 7r — /3 — 7 if one allows for the presence of 
NP in mixing. By now it is well known that the so-called methods to 
measure a in fact ought to be called methods 2 to measure (3 + 7 [23 12H] • Of 

2 Note that by definition a + (3 + 7 = tt and also a + /3 + 7 = 7T. This relation does 
not test 3x3 unitarity, as it cannot be violated since it is a definition. Recall that with 3 
complex numbers (the sides of the unitarity triangle - even if they do not close -) we can 
define only two independent relative angles: the third one is always a combination of the 
other two. 
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course, the representative example of this method as far as flavour content 
is concerned is B — > ttti. 

It has to be stressed that there are other available methods to extract 7, 
i.e. using b — > s transitions. However, in the search for NP, great attention 
has to be payed to the assumptions made in order to keep the errors under 
control. In this respect, SU(2) is a very good symmetry of the strong in- 
teractions, but the same does not hold for SU(3). Furthermore, there are 
no first principle calculations of SU(3) breaking effects in non-leptonic B- 
decays. Therefore, although analyses based on SU(3), [/-spin, can be quite 
useful in some instances, we will restrict ourselves to »S77(2)-based analyses, 
as far as strong interactions are concerned. 

Concerning the two previous methods, their fundamental differences lie in 
the presence or not of penguin diagrams. The method to extract 7 or 2/3 + 7 
relies on pure tree-level processes, therefore no additional assumptions are 
needed to get a value of 7 7^ so as to have an irrefutable proof for a complex 
Vckm- Of course it has to be understood that 2/3 has been extracted from 
B° d ,B°^ J/^K S , J/^K *, J/^K°* in order to use the method of the 2/3 + 7 
extraction. In the a method, the presence of penguins in general calls for the 
additional assumption of no NP competing with these penguins in the decay 
amplitudes. Therefore we will present two separate analyses with the actual 
relevant data. 

5.1 7 from pure weak tree-level decays 

The enormous effort developed at the B-factories Belle and BaBar has re- 
sulted in the first measurements of 7 - although by now still poor - in tree- 
level decays B ± -> DK^, B ± -> D*K ± -> (Dn ) K*, where the two paths to 
D° or D° interfere in the common decay channel D°,D° — > Kstc + tt~- From 
a Dalitz-plot analysis [23, Belle has presented 7 = 68° ± \f a ± 13° ±11° and 
BaBar [3D] 7 = 70° ± 26° ± 10° ± 10°, together with the solutions obtained 
by changing 7 — > 7 ± ir. The statistical significance of CP violation in the 
Belle measurement is of 98% [29J; by now therefore this can be taken as an 
indicative figure of the statistical significance of having a complex CKM. 

The method used by Belle and BaBar suffers from model-dependent as- 
sumptions in the Dalitz-plot analysis; these can be eliminated with more 
statistics following the suggestion of Ref. [21]. Also, the presence of NP 
in the D°-D° mixing could modify this analysis, but if necessary, in case 
of better precision, this small correction 24^ could be included in the way 
suggested in Refs. |3*Tj . 

We average conservatively both measurements to the value 7 = 69° ±21° 
(—111° ± 21°), which we take as a quantitative measurement of a complex 
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CKM matrix independent of the presence of NP at the one-loop weak level. To 
analyse the implications for the presence of NP, we add this new constraint 
to the previous analysis presented in Figs. |5(a)| and |5(b)| In Fig. |6(a)| we 
represent the analogue of Fig. |5(a)[ From Table ^ where we only present 
the central values, we clearly see that with two 7 values and two signs for 
cos(2/5) we generate four solutions. 




Figure 6: 68% (black), 90% (dark grey) and 95% (grey) probability regions, 
including the constraint on 7. 

It has to be stressed that although (3 presents a fourfold ambiguity, 2(3 
- necessary to extract 2<p d - presents a twofold one. With a measurement 
of 7 (another twofold ambiguity) we have j3 and R t fixed up to a twofold 
ambiguity, and therefore in total we have four different solutions for r^, 2(f) d 
and the CKM parameters. 

These solutions are represented in the (r^, 2<f) d ) plane as the joint PDF in 
Fig. |6(b)[ We have the "SM solution" corresponding to 2<p d ~ 0°, r d ~ 1, 
two NP solutions overlapping at 2<ft d ~ —85° and another NP solution near 
24>d ~ —169°. Had we used the constraint on the sign of cos(2/?), from 
B® — > J/^/K*, we would have eliminated two NP solutions, but the situation 
is not conclusive (see [H2])- The fact that the central value of 7 is a little bit 
high with respect to the value obtainable in a SM fit, implies a large value for 
Rt for the "SM solution", and therefore a low central value of r d , i.e. r\<\. 
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7 ( in P ut ) 


P 


20d 


T l 


Rt 


Sign[cos(2/?)] 


2/3 (input) 


a 


69° 


23.2° 


-0.8° 


0.80 


0.93 


+ 


47.2° 


87.4° 


69° 


23.2° 


-86.4° 


0.80 


0.93 




132.8° 


44.6° 


-111° 


-17.9° 


-83.0° 


0.48 


1.20 


+ 


47.2° 


-92.6° 


-111° 


-17.9° 


-168.6° 


0.48 


1.20 




132.8° 


-135.4° 



Table 1: "Central" values for the solutions (see Fig. |""J). 



Information on 2/5 + 7 is a l so available from partially reconstructed B° — > 
D*^ decays, |sin(2/3 + 7)| > 0.58 (90% C.L.) But at the moment 

this result relies on theoretical inputs based on SU(3), and we will therefore 
not use this constraint. Coming back to our irrefutable proof of a complex 
CKM, with our PDF for 7 we find that a real CKM matrix is excluded at 
a 99.84% C.L. 3 . In much more intuitive terms, we find that the probability 
corresponding to 7 G [10°; 170°] U [-170°; -10°] is 99.7%. 



5.2 7 from a. methods 

Using these a methods to get an irrefutable proof of a complex Vckm is m 
principle much more subtle. They are all based on the transition b — > duu, 
so that the so-called penguin pollution (weak loops) can also become NP 
pollution. The most prominent channels are B — > 7T7t, B — > pn and B — > pp. 
In the appendix we have shown that the extraction of a = vr — (3 — 7 = a + <pd 
from these channels is valid even in the presence of any NP in the AI = 1/2 
piece of the weak hamiltonian. As the AI = 3/2 piece is dominated by tree 
diagrams, this extraction of a is in accordance with our initial assumptions: 
"no NP in weak processes where the SM contributes at tree level" means in 
this case no NP in the AI = 3/2 piece. We conclude that the extraction of a 
from these channels is completely independent of the presence of NP in the 
usual penguins. Only NP in electroweak penguins (EWP) could spoil this 
conclusion, but this is very unlikely due to the small effect of the SM EWP, 
of order 1.5° (see reference [HE]), that can be eventually taken into account. 

BaBar has presented a time-dependent analysis of the p + p~ channel [3*6] . 
that once supplemented with the p + p° and p°p° branching ratios [27| EE] and 
the measurement of the final polarization [20], can be translated [SH] into 

3 With the standard construction used by Babar and Belle this C.L. is obtained as the 
integrated probability over the values of 7 having more probability than the most probable 
CP-conserving value of 7 |rS4| . 
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the measured value a = 96° ± 10° ± 5° ± 11° where the last error comes 
from the usual 577(2) isospin bounds 4 |a e g — a\ < 11°. Because the mea- 
surement is sensitive to sm(2a e ff), a = a e ff ± 11° presents a fourfold am- 
biguity (a, a + 7r, | — a, —a — | ) . In the pn channel the pentagon isospin 
analysis - from quasi-two-body decays - needs more statistics and/or ad- 
ditional assumptions. A time- dependent Dalitz-plot analysis in the chan- 
nel B — > 7T + 7r~7r° has been presented by BaBar [371 HO], with the result 
a = 113° ± f 7 °o ±6°. 

Since this analysis is sensitive to both sin (2a e ff) and cos(2a fr), the re- 
sulting ambiguity is just a twofold one (a, a + 7r). It is remarkable that these 
two solutions are in good agreement with two of the solutions coming from 
the pp channel. This important property will be used to eliminate two of the 
four solutions coming from the pp channel, in order to see the future trend 
of this analysis in a much clearer way. 

The situation in the nir channel does not yet allow a full isospin anal- 
ysis and the isospin bounds are quite poor. Furthermore BaBar and Belle 
measurements are still in some conflict, so that we will not use these results. 

As before, we average the data from pp and pn but only keep the two 
solutions consistent with the pit channel data. Our averaged values are a = 
100° ± 16°, (-80° ± 16°). 

The analogue of Fig. |6(a)| for the apex of the unitarity triangle, with the 
information on 7 replaced by that on a, is represented in Fig. |7(a)| 

In this case the data are not conclusive enough to eliminate a real Vckm] a 
real CKM is excluded at a 75% C.L.. In the second approach, the probability 
of 7 G [10°; 170°] U [-170°; -10°] is 85%. In Tabled where we include the 
central values, we have four solutions from the two a values and the two 
signs of cos(2/?). Again, the fourfold (3 ambiguity is effectively reduced to 
two because to reconstruct 7, the quantity that matters is (a + (3) =71 — 7. 
These solutions are represented in Fig. |7(b)| As before, we have the SM 
solution and the other three NP solutions that tend to overlap here more 
than in the previous case. 



Comparing the two tables and taking errors into account, it is easy to rec- 
ognize that the two solutions with positive sign of cos(2/3) can be consistent, 
but the negative ones are somewhat inconsistent. The reason can be seen by 
looking at the a generated from the 7 measurement. Because a = it — (3 — 7, 
the two a generated with the same (3 (the same sign of cos(2/3)) differ by 7r, 

4 In our notation cE e ff is the usual a e g of Refs. |25U14ll77j . but where we have introduced 
/3 instead of (3 as the phase in the B^-B^ mixing. 
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Figure 7: 68% (black), 90% (dark grey) and 95% (grey) probability regions, 
including the constraint on 7 through a. 



7 




2<j) d 


r l 


Rt 


Sign[cos(2/9)] 


2/3 (Input) 


OL (Input) 


56.4° 


22.8° 


-1.7° 


0.96 


0.85 


+ 


47.2° 


100° 


13.6° 


8.5° 


-115.7° 


1.78 


0.62 




132.8° 


100° 


-123.6° 


-15.1° 


-77.4° 


0.44 


1.26 


+ 


47.2° 


-80° 


-166.4° 


-3.8° 


-140.5° 


0.36 


1.39 




132.8° 


-80° 



Table 2: "Central" values for the solutions (see Fig. 0). 



so if one of the 7 measurements is consistent with one of the a measurements, 
so is the other, because both measurements of a and 7 differ by n. On the 
contrary, for a fixed 7 measurement, the two a generated with the different 
(3 (different sign of cos(2/5)) differ by 42.8°, so if one solution is compatible 
with the a measurement, the other is not. In our scenario, therefore, a way of 
measuring the sign of cos(2/3) relies on the measurements of sin(2/3), sin(2o7), 
cos(2a) and tan 7. 

In Fig. |H1 we repeat the previous analyses including both constraints 
from 7 and a. As can be seen the solutions with negative cos(2/3) have 
almost disappeared, in fact just 11.1% of the points have negative cos(2/5), 
contrary to the previous plots, where they were 50%. Combining 7 and a 
measurements, a real CKM matrix is excluded at a 99.92% C.L. and the 
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probability of 7 G [10°; 170°] U [-170°; -10°] is 99.7%. 




Figure 8: 68% (black), 90% (dark grey) and 95% (grey) probability regions, 
including the constraints on 7 derived from both direct 7 measurements as 
well as from a. 

As we have mentioned, previous work on these analyses of New Physics 
has been presented by Z. Ligeti in ^Hj- For similar experimental inputs and 
theoretical assumptions the results agree. 



5.3 Further comments 

The semileptonic asymmetry is A^l = Im (r 12 /M 12 ), where in our scenario 
M 12 is polluted by NP, M 12 = ±AM Bd e i2 ^, but - to the present level of 
precision - r i2 oc (a + be~ n + ce _l27 ) is not; a, b, c depend on QCD parameters 
and moduli of the first two rows of the CKM matrix [HJI321E]- Therefore, 
by including the measurements of AM^, aj/^K s and a determination of the 
sign of cos(2/?), Asl could be used in the future to obtain information on 
7. The impact of the actual value A SL = (-3 ± 7) x 10 -3 |Sj in Figs. El El 
and|Hlis very mild, but a future determination of Asl at the level of the SM 
precision (A^ = (-5.0 ± 1.1) x 10~ 4 |H] ) would be sufficient by itself to 
completely eliminate some 7 solutions. 

In some parts of our analyses we have used 3x3 unitarity. It is important 
to clarify that a direct measurement of 7 from tree-level decays (or in the 
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future from Asi) provides an irrefutable proof of a complex CKM matrix, 
independent of violations of 3x3 unitarity; this includes both 7 and a meth- 
ods. This is not the case of our analyses of NP, where we have extensively 
used 3x3 unitarity, for example in the extraction of (3 and Rt from moduli 
of the first two rows and 7 so as to obtain 20^ and r\. 

6 Summary and Conclusions 

We have carefully examined the present experimental evidence in favour of 
a complex CKM matrix, even allowing for NP contributions to ex, AM^, 
AM Bs , aj/^K S i B — > pp and B — > pir. First, we showed that, contrary to 
some statements in the literature [19], it is not feasible in practice to derive 
a non- vanishing value of \J\ using only the present information on moduli of 
the first two rows of the CKM matrix. We showed this by pointing out that 
a legitimate determination of \J\ from these moduli, although possible in 
theory, would in practice require a completely unrealistic high-precision de- 
termination of \V US \, \Vcd\- We then introduced a convenient parametrization 
of NP and examined the impact of AM# d and aj/^K s m obtaining experi- 
mental evidence of a complex CKM matrix. We then emphasized that the 
best evidence for CKM to be complex, arises from 7, either through a direct 
measurement or through a measurement of a, together with (3. We conclude 
that if NP does not pollute SM amplitudes dominated by tree level diagrams, 
a real CKM is excluded at a 99.92% C.L., and the probability of having 7 in 
the region [10°; 170°] U [-170°; -10°] is 99.7%. We also illustrated how the 
above measurements can be used to place limits on the size of NP, which is 
allowed by the present data. As we emphasized in the Introduction, having 
an irrefutable piece of evidence for a complex CKM matrix, in a framework 
where the presence of NP is allowed, has profound implications for models 
of CP violation. In the particular case of models with spontaneous CP vi- 
olation, a complex CKM matrix favours the class of models [H] where, 
although Yukawa couplings are real, the vacuum phase responsible for spon- 
taneous CP violation also generates CP violation in charged-current weak 
interactions. Conversely, the evidence for a complex CKM matrix, even al- 
lowing for the presence of NP, excludes the class of models with spontaneous 
CP violation and a real CKM matrix at a 99.92% C.L.. 
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\v ud \ 


0.9738 ±0.0005 


\v us \ 


0.2200 ± 0.0026 




0.224 ±0.012 


\Vcs\ 


0.97 ±0.11 


\v ub \ 


(3.67±0.47)10" 3 


\Vcb\ 


(4.13±0.15)10" 2 


O'J/^K.s 


0.734 ±0.054 


A SL 


(-3±7)10" 3 


7 


(69 ±21)° mod 180° 


a 


(100 ± 16)° mod 180° 




0.22000 ±0.00002 






l Vcd U3(a) 


0.21985 ±0.00002 


l Vcrf U3(b) 


0.21955 ± 0.00002 




(0.502 ±0.007)ps- 1 


m Bd 


5.2794 GeV 


fs d 


(0.20 ±0.03) GeV 


B B d 


(1.30 ±0.18) 


Vb 


0.55 


S (x t ) 


2.5745 



Table 3: Numerical inputs for the different calculations [2"T| I44j . 



Appendix 

In this appendix we explain why the extraction of a based in the so-called 
isospin analysis is valid even in the presence of New Physics in the AI = 1/2 
piece of the b — > d hamiltonian. This is true for B — > 7r7r, B — > pp and 
B — > pit. 

In the B — > titt case, as has been rephrased in [12], the basic ingredients 
to extract a are: 
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The full hamiltonian only contains A J = 1/2,3/2 pieces, and isospin 
is a good symmetry of final state interactions. With these ingredients 
one has 

X„ s (11) 

V ^+0 ^3/2 

and A + o is a "physical observable" that can be reconstructed (up to 
discrete ambiguities) with the directly measurable observables A + _ = 

-4^ and the branching ratios of the 6 available channels B l+ i — ► 7rV- J 
and their CP-conjugates; obviously A +0 = A(B + — > 7r + 7r ) oc A 3 / 2 and 

A+_ = A(B° -> 7r + 7r). 



(ii) If the AI = 3/2 piece of the hamiltonian is exactly the SM one (the 
tree level piece 5 ) we have 

A +0 = e" i2 V i27 = e +im . (12) 

It is important to stress that the extraction of A + o is valid in any model of 
NP that fulfills assumption (i), therefore the extraction of a from equation 
(j!2)l is valid in any model where the NP in the decay amplitudes only appears 
in the AI =1/2 piece. If full isospin analysis is not done, because the usual 
bounds are obtained assuming (i), the result is also valid. The B —> pp case 
is similar to the B — ► titt case. 

In the Dalitz plot analysis of B — > pir, the moduli and the relative phases 
of A+- = A(B° -> p + 7r-), A-+ = A(B° -> p-vr+) and A 00 = A(B° -> 
p°7T ), are measured, together with the CP-conjugate channels, and a global 
relative phase weighted by g/p. With a general weak hamiltonian H w with 
AI = 1/2, 3/2 pieces one obtains from isospin 

r -^3A2 _ - ^ + "+^ + +2A°° , , 

A 3/2 , 2 A+- + A-+ + 2A00' Udj 

where ^3/2,2 is the reduced matrix element of H W (AI = 3/2) with a I — 2 
final state. Assuming that H W (AI = 3/2) is the SM one, we have 

c -»2ff ^3/2,2 _ e -i2(0 +1 ) _ e i2a _ 
^3/2,2 

So the extraction of a from the Dalitz plot analysis in B — > pre is completely 
valid in the presence of any NP in the AI = 1/2 hamiltonian. 



inclusion of EWP gives a shift of 1.5° in a that we will neglect, although it can be 
included. 
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